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Spin triplet superconductivity with line nodes in Sr2Ru04 
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Several possible odd-parity states are listed up group-theoretically and examined in light of 
recent experiments on Sr2Ru04. Those include some of the /-wave pairing states, d(k) oc 
zk x k v (k x + iky) and z(k x — ky)(k x + iky) and other z(k x + ik y ) cosck z (c is the c-axis lattice 
constant) as most plausible candidates. These are time-reversal symmetry broken states and 
have line nodes running either vertically (the former two) or horizontally (the latter) , consistent 
with experiments. Characterizations of these states and other possibilities are given. 
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The possibility of the spin triplet pairing in Sr 2 Ru0 4 
is pointed-putEl soon after the discovery of the supercon- 
ductivity rP Experimental evidence for the spin triplet 
pairing is accumulated. The temperature-independent 
Knight shift for the external magnetic field parallel to 
the conducting plane (the basal plane) is thought to be 
the evidence for the sain triplet pairing with the d- vector 
aligned to the z-axisBP The fiSR experiment^ indicates 
that the magnetic field is spontaneously induced and the 
time reversal symmetry seems to be broken in the su- 
perconducting state. These experiments were explained 
by tJ^ p-wave spin triplet state with the order parame- 
terl 



d(k) = A z(k x + ik y ), 



(1) 



where An is a constant. This state is analogous to the 
Anderson-Brinkman-Morel (ABM) state, which is iden- 
tified as the A phase in superfluid 3 HeJ3^ The AMB state 
in superfluid 3 He has point nodes of the energy gap at 
k x = k y = on the spherical Fermi surface. How- 
ever, since the Fermi surface of Si'2Ru04 consists of three 
cylindrical surfaces, the above p-wave state has a finite 
energy gap on the Fermi surface. The previous experi- 
ments of specific heatH-' C(T) and NMR relaxation ratetf 
T{~ showed that about a half of the density of states re- 
mains at T — in the superconducting state. In order 
to explain this residual density of states, the non-unitary 
states ancLjthe orbital dependent superconductivity are 
proposed Ej'Eir 

However, it is reported in very recent experiments that 
the residual density of states are very small in the better 
samples- Le. C(T) cx T 2 and !Tf 1 cx T 3 at low temper- 
atures.li3Ef These temperature dependences are inter- 
preted as a consequence of the line nodes of the energy 
gap. The residual density of states observed in the pre- 
vious experiments may be due to the impurity effects on 
the anisotropic superconductivity with line nodes. 

An anisotropic energy gap model caused by the fi- 
nite range interaction in the two-dimensional plane is 
proposed to explain the line-node-like temperature de- 



pendences— pf the specific heat and relaxation rate in 
Sr9Ru04.li-$ The order parameter in this case is given 
by 



d(k) = Aoz(smak x + isinafcj,), 



(2) 



where a is the lattice constant in the basal plane of 
the tetragonal crystal D 4 ^. This anisotropic energy gap 
should give the exponential temperature dependence as 
exp(— oA/T) at temperatures below the smallest energy 
gap A in clean samples, where a is a constant of the 
order 1. As far as we know, however, there are no ex- 
perimental evidence for the existence of the small gap in 
Sr2RuC>4. n 

The previously classified pairing states^ with a line 
node either (A) break the four-fold symmetry in the 
basal plane such as xk x belonging to the two-dimensional 
representation E„, which is incompatible with the four- 
fold symmetric H C 2 behaviorjli 1 or (B) have the two 
component d vector such as a non-unitary bipolar sate; 
xk x + iykv which is incompatible with the Knight shiftcP 
and ^SRiS experiments. Thus, if we accept these ex- 
periments, all the previous states becomes unlikely to 
cxplainO This, is also true for those states listed by 
Rice and SigristcP and Sigrist and Zhitomirsky.EP 

Here we are going to first list up several possible odd- 
parity states with lower angular momentum allowed un- 
der the tetragonal symmetry D^h of Sr2Ru04, then we 
argue plausibility of these states in light of the existing 
data where a pairing mechanism is also examined for 
stabilizing some of the classified states. 

In Table I the orbital functions allowed in D^hi which 
are derived from the product of the two irreducible 
representations, are shown. These belong to the two- 
dimensional representation E„. These are used when 
constructing the order parameters, which are listed up 
in Table II. 



Among these allowed states the following /-pairing 
states are particularly attractive: 
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Table I. Orbital function in two-dimensional representation E v 
and its product representation 



Product 


Orbital function: (Ai(fc), X 2 {k)) 


A lg xE u 


(fe^ + ky)(k x , ky) 


A 2g xE u 


k x k y (k x — ky)(k x ,k y ) 


Bi 9 xE„ 


(k x — ky)(k x , k y ) 


B 29 xE„ 


k x ky (k x , ky ) 



(Aigiil.) k x -\- k^ A_2g'- kxky(k x ky), Bi^: k x T$2g- kxky~) 



Table II. Order parameter and its characterization 



order parameter 


Unitarity 


TRSB 


zAi(fc) 


U 


N 


(x + iy)Xi(k) 


N 


Y 


{x + iy)(Xi(k) + X 2 (k)) 


N 


Y 


2(Ai(fc) + A 2 (fc)) 


U 


N 


{x + iy)(Xi[k)+iX 2 {k)) 


N 


Y 


fcAi(fc) + yA 2 (fc) 


U 


N 


iAi(fc) + iyX 2 (k) 


N 


Y 



Unitarity (U) or Non-unitary (N), TRSB (time reversal symmetry 
breaking) yes (Y) or no (N) 



d(k) = A Q z(k 2 x 



k y )(k x "t" iky) 



and 



d(k) = Aozk x k y (k x + ik y ) 



(3) 



(4) 



because (1) they are both odd-parity states and have 
single d-vector component, and (2) the four vertical line 
nodes run along the c-axis and the four-fold symmetry 
in the basal plane is preserved. 

Here we should point out a possibility that these /- 
pairing states can mix p-wave pairing in general for D^, 
namely, 

d(k) = A z{k x k y (k x + iky) + ~/{k v + ik x )} (5) 

with 7 being a complex number. This mixing washes 
out the desired line node except for some special case. 
Within the present framework we assume that the an- 
gular decomposition is nearly complete, neglecting their 
mixing. At the moment we do not know how well this 
decomposition is, but the experimental facts seem to de- 
mand it. 

So far we only consider the purely two-dimensional or- 
der parameters characterized by the two wave vectors 
k x and k y in the basal plane. However, there is an- 
other class of the states compatible with the experimen- 
tal data, namely, 

d(k) = z[f (k x , k y )(ao + <X\ cos(cfc z ) + a 2 cos(2cfc z ) + ...) 

+fe(kx,ky)(bism(ck z ) + b 2 sm(2ck z ) + ...)], (6) 
where the order parameter in this class is not the partial 



wave in k z but the Fourier series with c being the lattice 
constant along the c-axis. fo(k x ,k y ) and f e (k x ,k y ) are 
the odd and even functions with respect to the inversion 
in the k x — k y plane, respectively. The previously consid- 
ered states (eq.([[|) and eq.(g)) are two-dimensional states 
with ao = 1 and a n — b n = (n > 0). We propose the 
state with a\ ^ and other a n and b n are zero. The 
order parameter is written as 



d(k) = Aoz(k x + iky) cos(cfc z ). 



(7) 



This particular state can be regarded as derived from 
a product representation Ai g xE u and is periodic along 
the c-axis. This lattice periodicity might be important 
because all the three Fermi surfaces are open along the 
c-axis, but close within the basal plane. In fact, this 
state has horizontal line nodes at k z = ±7?-, which run 
parallel to the basal plane, and breaks the time reversal 
symmetry. Even if the ao term is mixed, i.e. 



d(k) — A z(k x + ik y )(cos(ck z ) + a ), 



(8) 



the horizontal line nodes exist as long as ao is real and 
\a \ < 1. 

If we take account of the tetragonal K2N1F4 type crys- 
tal structure or the body centered tetragonal lattice of 
Ru, the order parameter, 

r/ "' (9) 



d(k) = Aoi:(sm + isin ■ 



-T-) COS — — 

2 2 ; 2 

is preferable. This state has the horizontal line nodes 
at k z = ±~. This state is possible if the effective in- 
tralayer interaction is repulsive and interlayer coupling 
is attractive, as we discuss below. 

Due to the Coulomb interaction the effective inter- 
action will be repulsive between electrons in the same 
plane. We assume the effective interaction is attractive 



for electrons at and v% ± | ± § ± |, 



V(r,r') 



-Vo, 



if r' = r ± % x ± § y ± 



(10) 



where Vo > 0, n and n ± | x ± |y ± |i are the lattice 
sites of Ru, and a and c are the lattice constant in the 
body centered tetragonal lattice. The Fourier transform 
of the effective interaction is 



V(k,k') = -SVbcos 
c(k 



a(k x 



x cos ■ 



2 

K) 



k' x ) a(k y - k') 
■ cos — 



. ak x 
-8Vo(cos — — cos • 



sin ■ 



akn 



■ sin ■ 



aki, 



. ak y nk ,,,,,, ■■"//■, 

x cos — - cos — - + sin — - sin — - ) 

y 2 2 2 2 ' 



x (cos — — cos — — 
v 2 2 



ck. 



ck~ 



sm sin 

2 2 



r) 



(11) 



The orbital part of the order parameter caused by this 
interaction should have the form 



sin ^f- cos ^j- cos : 2 
cos sm cos ■ 2 

aka, „ n „ afe a ck z 



ck z 
2 
ck z 



COS • 

sin : 



■ cos • 



2 vjvu 2 sm ■ 2 



(12) 
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in the triplet pairing. These four states have different 
transition temperature in the weak coupling limit, but 
the first two states are degenerate. Although we have 
to calculate in more detail to find which state is most 
stable, we think the order parameter 



d(k) = iAo(sin cos — 

ck. 



(13) 



will possibly be the stable state. If we neglect the k, x and 



k y dependences in cos and cos , we get the order 



parameter given m eq. 

We propose some experiments to test the horizontal 
and vertical line nodes. The thermal conductivity will 
have weak dependence on the direction in the x — y plane 
and have a power low dependence on temperature, if the 
line nodes are horizontal. On the other hand it will show 
strong four-fold symmetry and exponential dependences 
in some direction if the superconducting state has the 
vertical nodes of the energy gap. The penetration depth 
will also show the orientation of the line node. This 
kind of the orientation dependent transport measure- 
ments, including the ultrasound attenuation was deci- 
sive in identifying the intricate gap structure of another 
triplet superconductor UPtsO 

The Josephson effect between conventional supercon- 
ductors and Sr2RuC>4 is observed i-ajrd shows the in- 
terestiag .temperature dependences Jla^ Some explana- 
tions^' Ef are given. Although the experimental study of 
the direction-dependence of the Josephson effect will be 
difficult because of the roughness of the interface, the ex- 
istence of the Ru lamellas and 3K phases,EiP we consider 
the ideal Josephson junction. Since all the states con- 
sidered in this paper have d \\ z, the Josephson current 
I || z between the triplet and the singlet superconduc- 
tors isJjoxbidden in the first order in the spin-orbit cou- 
pling.E3E3 The Josephson current between spin triplet 
and spin singlet is possible, if the d vector is not parallel 
to the current and the total (spin plus angular) momen- 
tum is conserved. In the states given in eq. ([?]), eq.(||) or 
eq.(13), the Josephson current in the x — y plane is can- 
celed, since the order parameter changes the sign in the 
fc z -direction. However, the Fermi surface is corrugated 
as observed by angle-dependent magnetoresistance oscil- 
lationc5Ea) and. the angle dependence of the de Haas van 
Alphen effect fS 1 and as a result the cancellation will not 
be perfect. In the presence of the mixing of the k z inde- 
pendent term (as in eq.(||)) the Josephson current is also 
finite. 

The horizontal line nodes will be most clearly seen 
by the neutron scattering. Since the quasi-particles are 
excited around k z — ±- in the superconducting state 
eq.(||) or eq.(13), the excitations between these quasi- 
particles, i.e. excitations with q z — ±— will become 
large. The neutron scattering experiment below T c will 
have peaks at q z = ±— , while it is almost independent 
of q z above T c . 

As a final remark, we point out the following additional 
possibilities: n 

(1) Suppose that we disregard the fiSK experiment^ 



Among the previous states® the non-unitary bipolar 
state, described by 



d(k) = A (xk x + iyky) 



(14) 



has line nodes for both up-spin and down-spin pair 
branches, thus there is no residual T-linear specific heat. 
This state does not exhibit the macroscopic sponta- 
neous moment averaged over the Fermi surfaces below 
T c . Since the effective spin-orbit coupling felt by the 
Cooper pairs may be small in Sr2Ru04,tP the d vector 
can rotate in the spin space. When the external field 
H is in the basal plane, the non-unitary bipolar state 
becomes 



d(k) = A Q (zk x + \jk y ) 



(15) 



with j being a vector lying in the basal plane which is 
rotatable so as to keep H^L j, being consistent with the 
Knight shift experiment .EP As another possibility, the 
time reversal symmetry conserved state becomes a can- 
didate such as zk x or zky which are degenerate in D4/J,. 
They form a domain structure in real system which may 
preserve the overall four-fold symmetry. Note that there 
is no four-fold symmetric states without the variable k z 
in our classified states.tr 

(2) Suppose further that we disregard the Knight shift 
experiment^ in addition to the /iSR experiment .B Then 
the following singlet pairings with line nodes are another 
candidates, including 



A(k) 



Aftk x hy{k x 



for A 2g 



A(k) = A (kl - hi) for B lg 



A(k) = A k x k y for B 



2g- 



(16) 



(17) 



(18) 



These singlet states might be stabilize dB by the ob- 
served antiferromagnetic fluctuationscS' whose wave vec- 

In conclusion, we list up and examine the several 
possible superconducting states in Sr2Ru04, including 
some of the /-wave pairing states, zk x k y (k x + ik y ) and 
z(k x — ky){k x + ik y ) and other z(k x + ik y ) cosck z . These 
states explain the experiments of /xSR, NMR and the 
specific heat. Experiments to distinguish these states 
are proposed. 

The authors would like to thank Y. Maeno, K. Ishida, 
A. Sumiyama, M. Ichioka and A.G. Lebed for helpful 
discussions. 
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